options(warn = 1)
# power of original Lampedusa study

outcome1 <- data62$factor
run1 <- data62$lamp

power1 <- rdpower(data=cbind(outcome1,run1), samph=cbind(30,30))

## Power at half SD Tau: 24.1%
## Power at quarter SD Tau: 9.6%


# power of smallest sample replication

outcome2 <- data82$factor
run2 <- data82$twelve

power2 <- rdpower(data=cbind(outcome2,run2), samph=cbind(30,30))

## Power at half SD Tau: 99.99%
## Power at quarter SD Tau: 86.8%